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Comments on “On Stored Energies and 
Radiation Q” 

Miloslav Capek, Member, IEEE, and Lukas Jelinek 


The commented paper m claims to provide a new ex¬ 
pression for an energy stored around a general radiator. The 
major purpose of this comment is to show that the claim is 
unjustified. Alongside with this issue, it is pointed out that 
some of the core formulas of m are not completely correct, 
and that their correct form has in fact been derived elsewhere, 
though for the purpose of evaluating the quality factor Qz and 
not the stored energies. 

The major outcome of the commented paper ||T| is claimed 
to be equations (29)-(31), which should represent the total 
stored energy Wm + We, the stored magnetic energy Wm, and 
the stored electric energy We (using the notation of IT]), and 
their connection with (2), (4), (5) of fT]. Equations (4), (5) of 
III are however not valid under the assumptions of IT], see 
m and 0. 

In order to be more specific, we will first recall here a well- 
established result HI, m 


It follows from (2) that (2) of m can be rewritten as 
W^ + We= lim (w^ + ITe - -/Tad) = 

r^oo \ 1/ / 

/ 2 \ 

= \l I U\\H\f + e\\E\f-2j^y^dndT, 

(3) 

in which v = IjyfJIe. Equation ([^ together with (4), (5) of 
m thus results in 



(4) 

which is not generally correct as it lacks the last term on the 
RHS of Q, which is generally non-zero jSj Eig. 10]. This 
omission has in fact already been mentioned in 15] text above 
Eq. (72)] and 0 text above Eq. (4)], but has been missed by 
the author of 0 despite of 0 being cited in 0. 

Alongside with the above argumentation, it is worth noting 
that the questionable association (result of (4), (5) of 0) 



( 1 ) 


where X is the input reactance of the considered radia¬ 
tor, which is fed by a frequency independenj^ current Iq, 
and which is surrounded by an homogeneous, isotropic and 
non-dispersive medium of permittivity e and permeability 
Equation Q further assumes r to be the radial coordinate, 
dn = sin 9 d6 dip being the solid angle differential, E, H 
being the electric and magnetic intensity generated by the an¬ 
tenna, respectively, and F (r) = limj._>oo {rE (r) exp (jfcr)} 
with k being the wave-number. The total radiated power can 
be written as a, 0 

Fred = j \\F\f dn. (2) 

n 


= W^ + We, (5) 

4 ouj 

can be tested directly. To this point, the solid line in Eigs.[2|^ 
shows that dX/duj can be negative for radiating systems, 
unlike to the stored electromagnetic energy. This well estab¬ 
lished theoretical and experimental fact (see for example 0) 
is omitted by the author of 0. 

The next commentary concerns the seemingly erroneous 
relation of (29)-(31) of 0 with (4), (5) of 0 , which can 
be written as 


3 |j^|2 ^^(Vandenbosch) 


Lorjv 

Stt 

ujr/iy 

Stt 


Vl U2 


dp (r 2 ) 1 sin (fci?) 


f f f T* / \ (^2) 1 sm(kE) 


Vl V 2 


(6) 


where 
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^The assumption of the frequency independent current at the input port is 
essential for the validity of ^ and of d. 0 later on. This assumption is 
not used in (J. However, as me system at hand is linear, this input current 
normalization can always be made. Any result claimed in (T] as valid for 
general input current thus must be valid also for the frequency independent 
input current. It is also important to stress that the frequency independence of 
the input current does not imply the frequency independence of the current 
density on the entire radiator. 


^/^(Vandenbosch) _ 


/ (^* (ri) ■ J (r2)^ cos^R) 

Vl V2 

+ J J (p* (’’i) P (’’a) - (ri) ■ J (r 2 )^ sin (kR) dV 2 dVi 

Vl V2 

(7) 


is an energy defined in jh] Eq. (63), (64)], p (r) is the charge 
density, and rj is the wave impedance. In fact, the correct 
version of ([^ has previously been derived in 0, 0 and reads 
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(using the notation of this comment) 


^ _ ^^(Vandenbosch) 


Lorju 

Stt 

ujr]u 

Stt 


Re 


Vi Vs 


LJ 


p* (»’i) 


dujp{r 2 )\ cos (kR) 


duj j 


R 


dFi 




Vl v^2 

( 8 ) 

where the input current /q is assumed to be frequency indepen¬ 
dent. Formula differs from ([^ in the last two terms. The 
Figs.[T]|^show which of the two formulas better represents the 
frequency change of the input reactance calculated as a ratio 
between voltage and current at the input port. The comparison 
is performed on an example of thin dipole similar to that of m 
and on an example of Yagi-Uda antenna treated in 0. In order 
to further support our claims, an additional and independent 
evaluation of dXjdu) can be found. This evaluation relies on a 
combination of Q and results of ||9l Eqs. (4), (5)] or similarly 
ifTOl Eqs. (25), (26), (28)], which reads 


1 dX 

4 duj 
rjv 

IGtt 
rjv 


|y 1 2 _ ■p;^(Vandenbosch) 

J (ri)p(r2)}0(ri,r2) dV^adEi 


Vl Va 


1 


+ J(r2)}^(ri,r2) dy2dVi 


Vl V2 

2V M J duj 

n 


■ F* dO, 


lo 


(9) 


in which 


G (r'i,r2) = 


(II 


r-2|P 


) ji {kR) 


R 


( 10 ) 


with ji (kR) being the spherical Bessel function of the 1st 
order. The results of (|^ are also depicted in Eig. [T] and Eig. 
The error in (|^ observed in Eigs.[T][^is however not visible 
in ID, the reason most probably being the assumption of an 
unrealistic current distribution with a purely reactive input 
impedance at the driving port. In contrast, Eigs. m consider 
a rigorous full-wave solution. 

As a final comment, it is important to mention that the 
results of 0, ID, for example (j^, are not mentioned in 
m, although they address the related problem of the source 
concept of the quality factor Qz- 
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Fig. 1. The normalized frequency derivative of the input reactance of a thin 
dipole made of a perfectly electrically conducting strip of width w = L/200 
and of infinitesimal thickness. The dipole is fed in its center by a delta gap 
source with the input current normalized to 1 A. The quantity R denotes the 
input resistance of the antenna. The calculation has been performed in the 
FEKO full-wave integral equation solver ED. The frequency derivative of 
the input reactance is calculated from the port impedance (solid line), from 
(red dashed line), from jsj (blue dashed line), and from (green dashed 
(me). The dashed orange lines highlight resonances, while the full orange 
lines highlight antiresonances. 
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Fig. 2. The normalized frequency derivative of the input reactance of a thin- 
strip Yagi-Uda antenna (see the inset). The width of the strip and its thickness, 
as well as the meaning of depicted curves, are the same as in Fig.^ 
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